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ABSTRACT 



A tale of two cores: Triggered massive star formation in the 

bright-rimmed cloud SFO 75 
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Context. Bright-rimmed clouds (BRCs) are isolated molecular clouds located on the edges of evolved HII regions. Star formation within the 
f"^ BRCs may have been triggered through the propagation of photoionisation-induced shocks driven by the expansion of the HII region. 

Aims. The main focus of this paper is to investigate the current level of star formation within one of these clouds and evaluate to what extent, 
if any, star formation may have been triggered. 

Methods. We present a detailed multi-wavelength study of the BRC SFO 75, including 1.3 cm and 1.2 mm continuum, and I3 CO and ammonia 
Q ' spectral line observations. To build up a comprehensive picture of the local environment we complement our observations with archival data 
£j \ from the 2MASS, GLIMPSE and IRAS surveys. 

■ Results. The 13 CO and 1.2 mm emission reveals the presence of a dense core located behind the bright rim of the cloud which is approximately 
coincident with that of the IRAS point source. From an analysis of the IRAS and 1.2 mm fluxes we derive a dust temperature of ~ 30 K, a 
^ luminosity of L bol = 1.6 x 10 4 L G and estimate the core mass to be ~ 570 M Q . The higher resolution ammonia observations resolve the 1.2 mm 
T T ~J ■ core into two distinct cores, one directly behind the cloud's rim (Core A) and the second located slightly farther back (Core B). These have 
K> ' masses of 8-15 M Q and 3.5-7 M Q for Core A and Core B respectively, which are significantly larger then their virial masses. Comparing the 
^ morphology of Core A with that of the photon-dominated region and ionised boundary layer leaves little doubt that it is being strongly affected 
by the ionisation front. 2MASS and GLIMPSE archive data which reveal a small cluster of three deeply embedded (A v ~ 20 mag) high- and 
intermediate-mass young stellar objects towards Core A leads us to conclude that the star formation found towards this core has been triggered. 
In stark contrast, Core B appears to have a much simpler, almost spherical, morphology. No stars are found towards Core B. We find evidence 
supporting the presence of shocked gas within the surface layers of the cloud which appears to extend to midway between the two ammonia 
cores. 

Conclusions. The scenario that emerges from our analysis is one where the two ammonia cores pre-date the arrival of the ionisation front. Since 
its arrival the over-pressure of the ionised gas at the surface of the cloud has driven shocks into the surface layers of the cloud. The propagation 
of these shocks through Core A have triggered the formation of a small cluster of massive stars, however, the shock front has not yet propagated 
deeply enough into the cloud to have affected the evolution of Core B. 

Key words. Radio continuum: stars - Stars: formation - Stars: early-type - Stars: pre-main sequence - ISM: clouds - ISM: individual object: 
SFO 75. 
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1. Introduction 

Bright-rimmed clouds (BRCs) are isolated molecular clouds 
located on the edges of evolved HII regions. Star formation 
within the BRCs may thus have been triggered through the 
propagation of photoionisation-induced shocks driven by the 
expansion of their HII regions. The photoionisation of a BRCs 
surface layers by UV photons from nearby OB stars leads to 



the formation of a layer of hot ionised gas. This is known as the 
ionised boundary layer (IBL), which surrounds the rim of the 
molecular cloud. Shocks are driven into the molecular gas if the 
IBL is over-pressured with respect to the molecular gas within 
the BRC, resulting in the compression of the cloud. This can 
lead to the formation of dense cores, which are then triggered to 
colla pse by the same (or a subsequent) shock front dElmegreen 
1 1992b . This propagating shock front may also serve to trigger 
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the collapse of pre-existing dense cores, leading to the creation 
of a new generation of stars. 

This mode of star formation is known as radiative-driven 
implosion (RDI) and may be responsi ble for the product ion of 
hundreds of stars in each HII region ( Qgura et alj|2002l) . per- 
haps even contributing up to ~ 15 % of the low-to -intermediate 
mass initial mass function (Sugitani et al. 2000|). In order to 



evalua te t his mode of triggered sta r formation ISugitani et al 



dl991l) and lSugitani & Ogurald 19941) identified 89 BRCs (com- 
monly referred to as the SFO catalogue) where star formation is 
likely to be taking place. These clouds were identified by cor- 
relating IRAS point sources - having colours consistent with 
embedded protostars - with clouds displayin g optically brigh t 
rims from the Sharpless HII region catalogue (ISharplessll 19591) 
and the ESO(R) Southern Hemisphere Atlas. 



We have reported in an earlier paper the results of a set 
of radio continuum observations made toward s all 45 SFO ob- 
jects located in the southern hemisphere (see iThompson et al.l 
2004 for details; hereafter Paper I). In this paper we identified 
18 BRCs that are associated with centimetre continuum emis- 
sion which is spatially and morphologically consistent with the 
presence of an IBL (i.e., radio continuum emission is elon- 
gated parallel to the optical emission and lies in front of the 
cloud with respect to the direction of ionisation). Follow-up 
molecular line and high resolution radio continuum observa- 
tions were carried out on a sample of four clouds from Paper I 
in order to evaluate whether the presence of the ionised gas was 
having a direct influence on the evolution of these clouds and 
the inferred star formation within them. Dense molecular cores 
were found to be located directly behind the io nised rims with 



respe ct to the direction of their ionising stars (lUrquhart et al 



2006; hereafter Paper II). Evaluating the pressure balance be- 
tween the ionised and molecular gas, we identified two clouds 
in approximate pressure balance and therefore likely to be in a 
post-shocked state, and two clouds that are in the early stages 
of radiative-driven collapse with shocks still being driven into 
them. 

From Paper II we were able to identify two clouds at an 
early stage of interaction with the expanding ionisation front 
and found tentative evidence for star formation within these 
clouds. In this, the first of two papers, we present a detailed 
multi-wavelength study of SFO 75, with SFO 76 being the 
subject of a subsequent paper (Urquhart et al. 2007, in prep.). 
These include radio and millimetre continuum, I3 CO and am- 
monia spectral line observations as well as 2MASS, GLIMPSE 
and IRAS archival data. The main focus is to investigate the 
current level of star formation within this cloud and evaluate to 
what extent, if any, the star formation may have been triggered. 
The structure of this paper is as follows: in Sect. 2 we sum- 
marise the main features of SFO 75 and its host HII region; 
in Sect. 3 we describe the observational procedures and data 
reduction processes; we present the observational results and 
analyses in Sect. 4; our main findings are discussed in Sect. 5 
with particular emphasis on whether star formation has been 
triggered; and in Sect. 6 we present a summary of our results 
and highlight our main findings. 



2. SFO 75 and its environment: H a and 8 /urn 

SFO 75 is situated on the edge of the HII region RCW 98 that is 
located at a heliocentric distance of 2. 8 kpc. The HII region is 
ionised by LSS 3423, an 09.5 IV star dYamaguchi et al.ll 1 999t) 
which lies at a projected distance of ~ 0.8 pc from the op- 
tically bright-rim of SFO 75. Embedded within this cloud is 
the IRAS point source 15519-5430, which is the most lumi- 
nous in the SF O catalogue with an far-I R luminosity (Lfir) of 
~ 3.4 x 10 4 L m dSugitani & Qguralll994l) . 

In Fig. Q] we present a three colour composite image of 
SFO 75 and the HII region within which it resides. This im- 
age h as been constructed from an SuperCOSMOS H„ image 
(blue; IParkeret al.ll2005l) . optical DSS-POSS II image (red) 
and the Spitzer-IR AC 8 um (green) imag e obtained from the 



GLIMPSE survey (Benjam in et al.l 120031) . The H a emission 



arises from the recombination of ionised hydrogen within the 
HII region. The 8 yum band is dominated by emission from 
two polycyclic aromatic hydrocarbon (PAH) features centred 
at 7.7 |/m and 8.6 /mi. PAHs are excited by UV photons that 
manage to penetrate some way into the molecular material sur- 
rounding the HII region, and thus their presence marks the in- 
terface between the ionised gas within the HII r egion and the 
surrounding molecular gas (Leger & P uget 1984), commonly 
referred to as a photo-dominated region (PDR). 

The 8 /im and H ff emission regions can clearly be seen to 
be anti-correlated with each other except at the position of the 
BRC. The 8 fim emission forms a continuous unbroken, ap- 
proximately circular, ring (green) which completely surrounds 
the ionised gas (blue), indicating that the HII region is ionisa- 
tion bounded. The roughly circular appearance of the HII re- 
gion would suggest that it has expanded into a fairly homoge- 
neous, spherically symmetric, environment. In addition to the 
circular ring of emission, the 8 /im emission reveals the pres- 
ence of a dust lane running from the north-east to the south- 
west at an angle of ~ 45°, bisecting the HII region giving it a 
"coffee bean" like appearance. 

SFO 75 stands out dramatically against the otherwise clas- 
sical looking HII region. In Fig. 1 the bright rim is seen as a 
thin white region that curves slightly away from the ionising 
star. The thin rim seen in the optical plate is consistent with 
the cloud and the illuminating star being located either in the 
same plane of the sky, or with the cloud being located slightly 
in the foreground. However, the cloud cannot be too much in 
the foreground as we would expect to see a region of extinction 
in the H a emission from the dense gas within the cloud, and 
this is not observed. We therefore conclude that SFO 75 pro- 
trudes into the HII region and lies approximately in the same 
plane as the ionising star, and thus its projected distance from 
the ionising star is likely to be close to its true distance from it. 



3. Observations and data reduction procedures 

3.1. ATCA centimetre continuum and ammonia 
observations 

Observations of the ammonia (1,1) inversion transition and the 
1.3 cm radio continuum were made on the 29th April 2004 
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Fig. 1. Three colour composite image of SFO 75 and its associated HII region, RCW 98. The position of the IRAS point source 
and ionising star are indicated by a black and white cross. This composite image has been produced using a combination of 
the DSS red band optical image, SuperCOSMOS H„ image and GLIMPSE 8 //m image which are coloured red, blue and green 
respectively. 



using the Australia Telescope Compact Array (ATCA), which 
is located at the Paul Wild Observatory, Narrabri, New South 
Wales, AustraliaQ The ATCA consists of 6 X 22 m antennas, 5 
of which lie on a 3 km east-west railway track with the sixth 
antenna located 3 km farther west. This allows the antennas 
to be positioned in several configurations with maximum and 
minimum baselines of 6 km and 30 m respectively. 

The observations were made at two different frequency 
bands centred at 23,569 and 23,695 MHz (NH 3 (1,1)) using 
bandwidths of 128 and 8 MHz respectively. The 8 MHz band- 
width was divided into 512 channels, giving a frequency reso- 
lution of 15.6 kHz and a velocity resolution of 0.2 km s _1 per 
channel. The observations were made over a twelve hour pe- 

1 The ATCA is funded by the Commonwealth of Australia for op- 
eration as a National Facility managed by CSIRO. 



riod providing good hour angle coverage. To correct for fluctu- 
ations in the phase and amplitude caused by atmospheric and 
instrumental effects the total integration time was split into 
blocks of ~ 15 minutes of on-source integration time sand- 
wiched between two minute observations of the phase calibra- 
tor 1613-586. For absolute calibration of the flux density Mars 
was observed once during the observations for approximately 
ten minutes. To calibrate the bandpass the bright point source 
1253-055 was also observed once during the observations. The 
main observational parameters are summarised in Table [TJ 

The calibration and reduction of th e data were perfo rmed 
using the MIRIAD reduction package dSault et alJI 1995b fol- 
lowing standard ATCA procedures. The task UVLIN was used 
to subtract the continuum emission in the 23,695 MHz band- 
pass from the ammonia (1,1) spectra in the uv plane. These 
data were then CLEANed using a robust weighting of 0.5 to ob- 
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Table 1. Observational parameters for the ATCA NH 3 and cm Table 2. Observational parameters for the Mopra 13 CO obser- 
observations. vations. 



Parameter 


NH 3 (1,1) 


1.3 cm continuum 


Parameter 




Rest frequency (MHz) 


23,695 


23,569 


Rest frequency (GHz) 


110.201 


Total bandwidth (MHz) . . . 


8 


128 


Total bandwidth (MHz) . . . 


64 


Primary beam (arcmin) .... 


2.5 


2.5 


Number of channels 


1024 


Number of channels 


512 


32 


Velocity resolution (km s~') 


0.17 


Channel resolution (kHz) . . 


15.6 




Beam size (") 


34 


Velocity resolution (km s _I ) 


0.2 




Date of observations 


7th July 2005 


Date of observations 


29 April 2004 


Total integration time (hr) . 


-1.5 


Total integration time (hr) . 




12 


Obs. positions 


RA(J2000) Dec.(J2000) 


Phase centre 


RA(J2000) 


Dec.(J2000) 


Map centre 


15* 55'" 49.3 s -54° 39' 13" 




15* 55 m 49 s 


-54° 39' 13" 


Offset position 


16* 02'" 17.6 s -55° 19' 12" 



Array configurations EW357 

Flux density calibrator .... Mars 

Bandpass calibrator 1253-055 

Phase calibrator 1613-586 



tain the same sensitivity as natural weighting, but with a much 
improved beam-shape and lower sidelobe contamination. The 
CLEANed data were restored to produce an image data cube 
of the ammonia (1,1) emission. An image of the 1.3 cm radio 
continuum emission was produced from the second frequency 
bandpass. 

3.2. Mopra 13 CO observations 

In Paper II we presented CO observations made in June 2003 
towards SFO 75 with the Mopra telescope. Mopra is a 22 me- 
tre telescope located near Coonabarabran, New South Wales, 
Australia^ The telescope is situated at an elevation of 866 me- 
tres above sea level, and at a latitude of 3 1 degrees south. 

Emission maps were constructed from a number of sin- 
gle pointing observations made over a grid of positions in the 
12 CO, 13 CO and C ls O (7=1-0) rotational lines. These observa- 
tions revealed the presence of an embedded core within the rim, 
and allowed us to estimate its physical parameters (e.g., exci- 
tation temperature, density, mass), however, these maps proved 
too small (2' x 2') to study the kinematics of the molecular gas. 
We therefore re-observed this cloud in July 2005 using the re- 
cently commissioned on-the-fly mapping mode to map the dis- 
tribution of 13 CO in a 5' x 5' region centred on the embedded 
core identified from our previous observations. 

The receiver is a cryogenically cooled (~ 4 K), low- 
noise, Superconductor-Insulator-Superconductor (SIS) junc- 
tion mixer with a frequency range between 85-1 16 GHz, cor- 
responding to a half-power beam-width of 36"-33" (Mopra 
Technical SummarjQ). The receiver can be tuned to either sin- 
gle or double side-band mode. The incoming signal is separated 
into two channels using a polarisation splitter, each of which 
can be tuned separately allowing two channels to be observed 
simultaneously. The receiver backend is a digital autocorrelator 

2 Mopra is operated by the Australia Telescope National Facility, 
CSIRO and the University of New South Wales. 

3 Available at http://www.narrabri.atnf.csiro.au/mopra/ 



capable of providing two simultaneous outputs with an instan- 
taneous bandwidth between 4-256 MHz. 

The observational setup for our 13 CO map is summarised 
in Table|2] For observations of SFO 75 both polarisations were 
tuned to the I3 CO frequency, with the the second channel be- 
ing periodically tuned to 86.2 GHz (SiO maser frequency) to 
check the telescope pointing at hourly intervals. The point- 
ing offsets were found to be < 10" rms. A total of 35 rows 
were scanned to produce the map with an off-source position 
being observed at the completion of each row to allow sub- 
traction of sky emission. Each row was scanned at a speed of 
~ 3.2 arcsec s _1 with the data being averaged over a 2 second 
cycle time. Individual rows were separated by 9". To correct 
the measured antenna temperatures (7^) for atmospheric ab- 
sorption, ohmic losses and rearward spillover, a measurement 
was made of an a mbient load (assumed to be at 290 K) follow- 
ing the method of Kutner & Ulich dl98lh approximately every 
15 minutes. System temperatures were found to be stable dur- 
ing the observations; r sys ~ 500 K with variations of no more 
than 10% during the map of SFO 75. Absolute calibration was 
performed by comparing measured line temperatures of Orion 
KL and M17SW to standard values. We estimate the combined 
calibration uncertainties to be no more than 20%. 

The data were reduced using LIVEDATA and GRIDZILLA 
packages available from the ATNF. LIVEDATA performs a 
bandpass calibration for each row using the off-source data 
followed by fitting a user-specified polynomial to the spec- 
tral baseline. GRIDZILLA grids the data to a user-specified 
weighting and beam parameter input. The data were weighted 
using the periodic r sys measurements and a cell size of 12" was 
used to grid the data. 

3.3. SEST millimetre continuum observations 

SFO 75 and its host HII region were observed in August 
2002 using the 37-channel bolometer array SIMBA (SEST 
Imaging Bolometer Array) at the Swedish ESO Submillimetre 
Telescope (SEST). SIMBA operates at a central frequency of 
250 GHz which corresponds to a wavelength of 1.2 mm, with 
a main beam efficiency of 0.5 and bandwidth of 50 GHz. The 
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beam has a half power beam width of ~ 24" for a single ele- 
ment with a sky separation of 44" between elements 

An area of approximately 800 arcsec 2 surrounding the 
IRAS point source associated with SFO 75 was mapped 
using the fast mapping mode of S IMBA (for details see 
Weferlin g et alJl2002l: iReichertz etai]|200lh . Two maps of 30 



minutes' duration and with a scan speed of 80 arcsec s were 
coadded to make the final map. The data were reduced follow- 
ing the standard procedures of the MOPSI reduction package, 
which was written by Robert Zylka (Grenoble Observatory j§. 
Atmospheric extinction was measured by skydips before and 
after each map, with zenith opacities at 225 GHz of typically 
0.3. Pointing was checked on an hourly basis and was found to 
be within 5". Absolute flux calibration was carried out by mea- 
surements of the 1.2 mm flux of Uranus. The final rms noise 
level in the map was found to be ~ 20 mJy beam 1 and the 
calibration uncertainty is estimated to be ~ 10%. 

4. Results and analysis 

4.1. Image analysis 

We present the images obtained from our observations in Fig. [2] 
The centimetre continuum emission can clearly be seen to form 
a thin layer that is elongated parallel to the optical rim and lo- 
cated on the outside edge of the PDR, thus confirming the pres- 
ence of an IBL between the HII region and the molecular cloud. 
The 13 CO emission traces the molecular gas which is concen- 
trated behind the bright rim falling off steeply in the direction 
of the HII region. The lowest 13 CO contour is almost parallel 
with the major axis of the IBL and the bulk of the molecular 
material is elongated parallel to the bright rim and perpendicu- 
lar to the direction of ionisation, consistent with the hypothesis 
that the cloud is being compressed by the ionisation front and 
its associated shocks. 

The 1.2 mm continuum emission appears to be elongated 
parallel to the ionisation front with a compact core of bright 
emission located almost directly behind the optical rim. This 
emission arises from a combination of optically thin thermal 
dust emission and thermal free-free emission from ionised gas. 
It is therefore possible that a significant proportion of the ob- 
served emission does not arise from the presence of an embed- 
ded core but from ionised gas associated with the bright rim. 
However, we can estimate the contribution from ionised gas at 
1.2 mm by extrapolating from the integrated 1.3 cm emission 
(30.3 mJy; see Paper II for details) assuming a thermal spectral 
index, i.e., v" ', which we find to be ~ 24 mJy. This is less than 
one per cent of the total integrated emission of 7.1 Jy. The con- 
tribution to the total flux from the thermal free-free emission is 
negligible and therefore the measured flux predominately arises 
from thermal dust emission. The 1.2 mm emission has a sim- 
ilar distribution to that of the 13 CO emission, although not as 
extended and peaks closer to the bright rim than the 13 CO emis- 
sion (~ 26"), possibly indicating that the SIMBA observations 



4 For more details see http://www.ls.eso.org/lasilla/Telescopes/SEST 
/html/telescope-instruments/simba/index.html. 

5 See also http : //www . astro . ruhr-uni -bochum . de /ni elbock/ 
simba/mopsi . pdf 



are sampling warmer material located closer to the rim. The 
13 CO and dust emission peaks appear to be centrally condensed 
and are located slightly south of the IRAS position. 

The higher angular resolution of the ammonia data resolves 
the 1.2 mm core into two distinct components, one that is lo- 
cated directly behind the bright rim (hereafter Core A), and a 
second which is located farther back from the rim (hereafter 
Core B) directly behind Core A with respect to the direction of 
ionisation. The face of Core A nearest to the bright rim has a 
flattened morphology which almost exactly matches that of the 
inside edge of the PDR and the radio emission. This morpho- 
logical correlation clearly indicates that these three regions are 
interacting with each other and that the morphology of Core 
A is being directly influenced by the UV illumination of LSS 
3423. Core B appears to have a more spherical morphology, 
possibly indicating that it has not so far been significantly in- 
fluenced by the passage of the photoionisation-induced shocks. 

In Table 3 we present a summary of the core parameters 
determined from the various emission maps. 

4.2. Analysis of global properties of the ionised and 
molecular gas 

In Paper II we presented a detailed analysis of the internal and 
external gas found towards SFO 75 using our previous CO and 
1.3 cm radio continuum observations, we therefore only pro- 
vide a summary of the relevant findings here. From our CO 
analysis we found that SFO 75 has a velocity of -37.5 km s _1 , 
a core kinetic temperature of ~ 24 K, a particle density of 
n = 3.4 x 10 4 cuT 3 and a core mass of Mb co = 177 M Q . 
Analysis of the ionised gas determined the peak ionising pho- 
ton flux to be 260 x 10 8 cirT 2 s _I and electron density to be n e 
~ 839 cm 4 . We estimated the internal and external pressures 
and found that the pressure of the ionised gas is approximately 
three times higher than the internal pressure of the embedded 
CO core. 

A drawback with our previous CO observations was the 
relatively small size of the maps (2' x 2') and rather coarse 
sampling of the data which limited our analysis, particularly 
with respect to investigating the kinematics of the cloud. For 
this reason we re-observed this cloud using the 13 CO (7=1-0) 
line, mapping a larger region (5' x 5')- The integrated emis- 
sion is shown in the right panel of Fig. [2] In Fig. [3] we present 
a position-velocity diagram taken through the cloud along the 
star-CO core vector (as indicated in Fig. [2]) starting at the po- 
sition of the ionising star (north-west corner of the image) and 
cutting through the core, running right to left through the dia- 
gram. 

This diagram reveals the molecular material to have a steep 
emission gradient on the surface exposed to the ionising star, 
and exhibits large velocity wings (~ 6 km s~') which are 
perpendicular to the direction of ionisation. These are indica- 
tive of shocked gas, with the blue shifted component being 
associated with compressed gas moving towards us and the 
red shifted component being associated with co mpressed gas 
moving away from us dLefloch & Lazarefnll994t) . The simi- 
larity between the position-velocity diagram predicted by the 
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Fig. 2. Left panel: Large scale 8 pan GLIMPSE image of SFO 75 and its associated HII region, RCW 98. The position of the IRAS point source and ionising star are indicated 
by a white cross. At 8 /urn the emission is dominated by PAH emission that are found in photo-dissociation region and marks the interface between the ionised gas within 
the HII region and the surrounding molecular gas. Right panels: Observational results laid over the region of the 8 fim GLIMPSE image outlined in white in the left panel. 
Clockwise from the top left panel: 1.3 cm radio continuum, integrated 13 CO emission, integrated ammonia (1,1) emission (main line) and the 1.2 mm continuum emission. 
The position of the IRAS point source is indicated by a white cross and its positional uncertainty by a black ellipse. The respective beam shape of each telescope is shown to 
size in the lower left corner of each plot. Additionally, in the image showing the 13 CO emission we have plotted a dashed line to indicate the cut through the cloud used to 
produce the position-velocity diagram presented in Fig. [3] The integrated 13 CO and ammonia emission has been integrated between -32 and -40 km s _1 , and -34.5 and -40 
km s _I velocity ranges respectively. The contour levels used in the image are as follows: 1.3 cm radio contours begin at 3cr and increase in steps of 2cr where cr = 2.5 mJy; 
13 CO contours begin at Act and increase in steps of 3cr where cr - 9 K km s ; 1.2 mm contour levels are 0.1, 0.2, 0.4, 0.7, 1.1 and 1.6 Jy; ammonia contours begin at Act and 
increase in steps of 3cr where cr = 80 mJy km s . 
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Table 3. Image parameters and peak positions of detected cores. 



Core tracer 


Beam parameters 


Position 


IRAS offset 


Deconvolved core 




ma j x min (") 


PAD 


RA(J2000) 


Dec.(J2000) 


(") 


FWHM (") 


"CO 


33 




15" 55"' 50.9 s 


-54° 39' 26" 


28 


41 


1.2-mm 


24 




15* 55 m 48.7 s 


-54° 39' 8.4" 


18 


36 


NH 3 (Core A) 


11.5 x 4.9 


-9.2 


15* 55 m 47.5 s 


-54° 39' 10.2" 


27 


10.4 


NH 3 (Core B) 


11.5 x 4.9 


-9.2 


15* 55 m 49.5 s 


-54° 39' 16.6" 


20 


6.6 




I I I I ^ 

-100" -50" 0" 50" 100" 

Angular Offset 

Fig. 3. Position-velocity diagram taken through the ionising 
star-CO core vector as indicated in Fig. [2] The direction of ion- 
isation is from right to left. This diagram clearly reveals the 
presence of blue and red shifted emission at the front of the 
molecular cloud, indicative of shocks and the steep emission 
gradient consistent with compression of the cloud by the ioni- 
sation front. The contours start at 20% and increase in steps of 
10% of the peak value (1 1. 1 K km s" 1 ). 




Log 10 (\[/im]) 

Fig. 4. Spectral energy distribution of SFO 75. Results of a two 
component greybody fit to the IRAS and 1 .2 mm dust contin- 
uum emission. 



similarities between our observational results and the predic- 
tions of the RDI model. 

The shocked gas layer can be clearly distinguished from 
the bulk motion of gas behind by its significantly broader ve- 
locity distribution. We are therefore able to estimate how far 
this shocked layer extends into the cloud from the FWHM of 
the shocked region. We thus estimate the shocked region to be 
~ 30", which corresponds to a physical distance of ~ 0.4 pc, 
extending from the edge of the cloud's rim to just in front of 
the CO core. Interestingly, this would place it approximately 
midway between the two ammonia cores as they are seen in 
projection. 



Lefloch & Lazarefn (119941) RDI model and our observations is 
striking (c.f. with their Fig. 14a). The presence of an IBL at 
the surface of the cloud, detection of a dense embedded core 
and the large pressure difference between the two regions led 
us to conclude in Paper II that this cloud was likely to be in 
the early stages of radiatively driven collapse. This conclusion 
is strengthened with the detection of the shocked gas and the 



4.3. 1.2 mm emission analysis 

4.3.1 . Luminosity, dust temperature and mass of the 
mm core 

In Fig. 2 we presented a contoured image of the 1.2 mm dust 
emission detected towards SFO 75, which has a peak flux den- 
sity of 2.0 Jy beam 1 and integrated flux of 7.1 Jy. In this sub- 
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section we use these flux measurements to derive the luminos- 
ity, temperature and mass of the embedded mm core. 
It has been shown (e.g., Sridharan et alj2002 ; Beuther et al 



2002) that two components are present in the spectral energy 
distribution (SED) of luminous YSOs: a compact hot compo- 
nent which dominates the IRAS 12 fim and 25 fim fluxes, and a 
more extended component arising from colder gas which dom- 
inates the IRAS 60 fim and 100 //m fluxes. We therefore es- 
timate the luminosity and temperature of these components by 
fitting two greybody function s to the core's SED . The greybody 
functions have the form (e.g. jDent et al.ll 1998b . 



F v = QB v (r dust )(l-e- T >), 



(1) 



where F v is the flux measured at a frequency v, Q is the solid 
angle subtended by the cloud, B v (Tdust) is the Planck function 
evaluated at the dust temperature (Td lls t) and frequency (v), and 
t v is the optical depth at frequency v. 

The SED of the 1.2 mm core was derived using fluxes from 
the four IRAS bands and the integrated 1.2 mm dust emis- 
sion from the SIMBA map shown in Fig. [2] The optical depth 
t is parameterized in terms of a common dust emissivity in- 
dex p so that t can be evaluated at arbitrary frequencies from 
a known reference frequency v re f and optical depth r re f, i.e., 
t v = T re f(v/v, e fy 3 . p has typical values of 1-2 for molecular 
clouds, however, to limit the number of free parameters, and 
to maintain consistency with previous wo rk of this kind (e.g., 
Faundez et al. 2004; Sridharan et al ] l2002h . we have adopted a 
value of 1 for the hot component, but leave the emissivity of 
the cold component as a free parameter. We also assume that 
the size of the region giving rise to the cold component emis- 
sion is the same as the 1.2 mm core. We leave the temperatures 
of the hot and cold components, and the size of the hot com- 
ponent as free parameters. The two component greybody fit to 
the SED is presented in Fig.Hl the temperatures for the hot and 
cold components are 175 K and 30 K respectively. We estimate 
the size of the hot central region to be ~ 0.2 " and the dust 
emissivity of the cold component to be ~ 1 .5. It is the cold dust 
temperature that is of most interest to us as it best describes the 
temperature of the large scale emission traced by our 1 .2 mm 
observations. 

Integrating over the observed SED yields a bolometric lu- 
minosity of ~ 1.6 x 10 4 L , which corresponds to the presence 
of a s ingle ZAMS s tar with a spectral class between B0 and 
B0.5 dPanagialll973h . Alternatively, if the observed emission is 
due to the presence of an embedded cluster, the luminosity of 
the cluster's most massive member would typically be one or 
two subclasse s lower than assum ing a single star, i.e., B2 or 
B3 star dWood & Churchwelll[l989l) . Since the majority of the 
luminosity is provided by the IRAS 60 fim and 100 fim flux 
densities it is necessary to consider where the emission origi- 
nates. Taking into account the size of the IRAS beam (~ 2' at 
100 fim), which is much larger than the mm core, it is possible 
that some of the measured flux arises from the ionisation front 
as well as the embedded core, and therefore the derived lumi- 
nosity should be considered as an upper limit. In order to better 
constrain the SED we attempted to obtain flux density mea- 



surements from the MIPSGAL Spitzer Legacy project which 
has recently finished surveying the GLIMPSE region at 24 fim 
and 70 /mi. However, this particular source lies on the edge of 
one of the observed fields, making it impossible to extract use- 
ful flux density measurements. Given the location of the IRAS 
point source and its close proximity to the mm core (~ 18") we 
would suggest that the majority of the flux is associated with 
the core, and that the core's true luminosity is close to the up- 
per limit. 

As stated previously, the majority of the millimetre emis- 
sion is due to optically thin thermal dust emission and by mak- 
ing a few assumptions it is possible to relate the integrated mil- 
limetre flux densit y to the total mass (gas+dust) of the core via 
dHildebrandll 19831) . 



Fl.lmmD 2 



M (gas+dust) = g R mm Mq, (2) 

^1.2mmOI.2mm(i dust/ 

and Ki.2mm are the integrated flux density and 
dust opacity per unit mass measured at 1.2 mm respectively, 
fii.2mm(7dust) is the Planck function for a dust temperature Tdust, 
D is the distance in kpc and g is the gas-to-dust ratio, which 
we assume to be 100. Using a value for /ci.2mm = 1 cm 2 g -1 
( Ossenkopf & Henningl ll994) and a value of 30 K for the dust 



temperature (as determined for the cold dust component of 
the greybody fit) we estimate the mass of the mm core to be 
~ 570 M . Given the uncertainties involved in this calculation, 
and that of the CO derived mass (i.e., Mu co = 177 M ), we 
find them to be in reasonable agreement. 

The luminosity and temperatures of the hot and cold com- 
ponents, which are similar t o those found toward s other mas- 
sive embedded objects (e.g., ISridharan et al]|2002l) . and the de- 
rived mass leads us to conclude that massive star formation 
is currently taking place within SFO 75. This is consistent 
with the statistical trend for BRCs to preferentially form more 
massive stars or multiple stellar systems t han found towards 
Bok globules and isolated de nse cores ("e.g.. lSugitani & O gura 
ll994ilYamaguchi et alJll999h . 

4.3.2. Molecular ring surrounding the HII region 

The only 1.2 mm emission detected in our wide-field SIMBA 
observations was that associated with SFO 75. This implies that 
the molecular shell of material surrounding the HII region (in- 
ferred from the presence of the PDR traced by the 8 fim PAH 
emission) must be of relatively low column density. We are 
therefore unable to determine the gas density within this ring, 
however, we can use the image rms to place an upper limit 
on the mass/beam. Assuming that the total mass is contained 
within a cylinder with a physical radius equal to that of the 
telescope beam on the sky and a length assumed to be equal to 
the width of the PDR seen in projection (i.e., ~ 1 pc) we can 
estimate an upper limit for the density of the surrounding gas 
to be 10 3 cirT 3 . 

This upper limit for the H2 particle density in the ring 
implies that the expanding HII region has not yet swept up 

6 For details see http://ssc.spitzer.caltech.edu/legacy/ 
abs/carey .html . 
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different hyperfine components. These profiles were fitted us- 
ing the CLASS software package, in particular the METHOD 
NH3 utility, assuming that all of the components have equal 
excitation temperatures and that the line widths and separa- 
tions are identical to laboratory values. Once the optical depth 
is known the excitation temperature of the ammo nia (1,1) in- 
versio n can be derived using Eqn. 1 and 2 from lHarju et alJ 
(1993) 




r B (l, 1, m) = r,^ (F(T eK ) - F(7b g )) (l - e- T(u ' m) ) , (3) 



-40 -30 
VLSR (km/s) 



Fig. 6. Spatially integrated spectral line profiles for Core A and 
Core B. The vertical dashed-dotted line indicated the velocity 
of the 13 CO core. 



where the function F(T) is defined as 
F(T) ' 



e hv/kT _ 1 ' ' 

Tb is the main beam brightness temperature, T ex is the source 
excitation temperature and 7b g is the cosmic background tem- 
perature assumed to be ~ 2.7 K, rj is the beam filling factor, 
which we assume to be of order unity, and m refers to the 
main hyperfine group. We have estimated the colum n density of 
mole cules in the upper (1,1) state using Eqn. 3 from lHarju et al.l 
dl993h . 



much of the surrounding material, either because the HII re- 
gion is relatively young or that the surrounding molecular gas 
is of low density. The low density of the molecular shell sur- 
roun ding RCW 98 suggests t he "collect and coll apse" process 
(e.g., De harveng et al.ll2005t IZavagno et al. 2006) has not oc- 
curred around RCW 98 and that SFO 75 is likely to be a pre- 
existing molecular cloud recently exposed to the UV illumina- 
tion of LSS 3423. 



4.4. Analysis of the ammonia cores 

In Fig. [5] we present eight velocity channel maps covering the 
velocity range between -33 to -38 km s . These channel 
maps highlight important differences between the two ammo- 
nia cores. First, the two cores are distinct both spatially and 
in velocity, with Core A and Core B centred at approximately 
-37 km s _1 and -35 km s _1 respectively (see also Fig. 6). 
Secondly, they illustrate the excellent correlation of the mor- 
phology of the flattened forward face of Core A with that of 
the PDR traced by the 8 /im emission - the morphologies of 
these two regions leaves little doubt that the ionisation front 
is having a dramatic influence on the evolution of this core. 
Additionally, the ends of the flattened face appear to curve into 
the cloud, reminiscent of the small "ears" seen to develop in 
the lLefloch & Lazareffl d 1 994 ) model (see their Fig. 4b) where 
the compressed layer meets the rear edge of the initial conden- 
sation. Conversely, Core B appears to have a more spherical 
morphology indicating that the shock front has not yet reached 
it. This is confirmed by the position of the shock front which is 
found between the two cores A and B (see Sect. 4.2). 

In Fig. [6] we present spectral line profiles integrated over 
each ammonia core. The hyperfine components can clearly be 
seen in these two spectra and can thus be used to estimate the 
optical depth for each core using the intensity ratio between the 



1.6 x 10 13 F(r ex )Av t(1,1,w) cm -2 , 



(5) 



where Av is the FWHM of the main line in km s _1 . To ob- 
tain the total column density for molecules in the (1,1) state we 
need to sum the populations of both the upper and lower states. 
This can be done by applying the Boltzmann equation assum- 
ing both the upper and lower levels are evenly populated, 



N(l, I) - N u + N\ = N u (l + e /,v/w ») . 



(6) 



The following equation calculates the total column density 
Af(NH3) making the assumption that only metastable levels are 
populated, 



-101.5/7-12 



N(NH 3 ) = N(l, l)^e 23 - 4/r ' 2 + 1 + ^ e - 4L5/ri2 + ye 

(7) 

Since there is no data available for the ammonia (2,2) tran- 
sition we are unable to calculate the rotational temperature, T12. 
However, we are able to calculate a range of values for each pa- 
rameter of interest, such as column and particle densities and 
core masses, using a range of rotational tempe ratures typical 
of similar cores (10^-0 K; e.g. JWu et al.ll2.Q06b . This range of 
temperatures corresponds to a difference of at most a factor of 
two in column density and so will not affect our results signifi- 
cantly. In order to obtain H2 column densities we have assume d 
an ammonia fractional abundance of 3 x 10~ 8 (IWu et al.ll2006l) . 
The H2 densities and core masses have been calculated using 
the core diameters presented in Table [3] The results of the fits 
to the ammonia spectral profiles and core-averaged values for 
the derived physical parameters of the cores are presented in 
Table [4] We consider the masses and densities of the two cores 
to be lower limits due to the missing short spacings data; tak- 
ing into account the uncertainties involved in their estimation 
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Fig. 5. Velocity channel maps of the ammonia (1,1) main line emission contoured over an enlarged region of 8 fim image 
presented in Fig. [2] (highlighted white square in lower right image). The central velocities are indicated in the top right corners 
of each map and are 0.6 km s" 1 thick. The position of the IRAS point source and its associated error are indicated by a cross and 
ellipse respectively. The first contour starts at 4<x and increase in steps of 4<x. The size of the synthesised beam is shown to scale 
in the lower left corner. 



Table 4. Fitted and derived parameters for the ammonia cores. 



Id 


Vel. 


T 

1 ex 


FWHM 


T 


N(H 2 ) 


«(H 2 ) 


Mass 


Virial Mass 




(km s" 1 ) 


(K) 


(kms- 1 ) 




(cm- 2 ) 


(cm" 3 ) 


(M G ) 


(Mo) 


A 


-37.0 


8.0 


1.63 


1.41 


4.0-7.7 x 10 22 


0.7-1.6 x 10 5 


8-15 


2.0 


B 


-35.3 


6.4 


1.23 


2.58 


4.3-8.6 x 10 22 


1.5-3.1 x 10 5 


3.5-7 


0.6 



they are likely to be significantly larger. Both cores appear to be 
centrally condensed with peak temperatures and densities a fac- 
tor of ~ 1 .5 and 2 larger, respectively, than the core-averaged 
values presented in Table [4] We can test the stability of these 
cores by comparing their masses to their virial masses, which 
we calculate using the standard equation (e.g., Evansfl999l) . 



=<Av> 2 , 



(8) 



where R C0K is the ammonia core radius in parsecs and Av is the 
FWHM line width. The resulting virial masses for each core 
are presented in Table [4] Comparing these to the derived core 
masses reveals that both cores are significantly more massive 
than could be supported by their internal pressure alone (in- 
ferred from their line widths). Therefore these cores are grav- 
itationally bound and given the large difference between their 
actual masses and their virial masses it is likely they are in a 
state of gravitational collapse. 

The visual extinction can be estimated using the column 
density (Nu 2 ) towards each core using the following relation- 



ship: A v = Nh 2 /0.94 X 10 21 mag dFrerkinget al.lll982l) . This 
results in similar values for both cores, A v ~ 42-82 and 45-91 
for Core A and Core B respectively. The visual extinction can 
be converted to K-band extinction (A^) by multiplying by the 
ratio between visual and infrared extinction i.e., A# s = A^/8.9 
( Rieke & Lebofsky 1985b . In this way we obtain K-band ex- 
tinctions between five and ten magnitudes for both cores. 



4.5. 2MASS and GLIMPSE point sources 

In this subsection we present complementary near- and mid-IR 
archival data obtained from the 2MASS and GLIMPSE surveys 
to investigate the stellar content of the ammonia cores. In the 
left panel of Fig. [7] we present a three-colour composite image 
produced by combining the 2MASS H, J and K s band images 
coloured blue, green and red respectivelyQ To show the spa- 
tial relationship between different regions (i.e., the dense cores, 
PDR and any embedded sources detected in the 2MASS data) 

7 Obtained from http : //www. ipac . caltech.edu/2mass/ . 
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Fig. 7. 2MASS and GLIMPSE three colour composite images. Left panel: 2MASS composite produced by combining J, H and 
K s band images which are shown in blue, green and red respectively. This image is overlaid with contours of the GLIMPSE 8 fim 
emission and the integrated NH3 emission in orange and blue respectively. Right panel: GLIMPSE composite made from the 3.6, 
4.5 and 5.8 //m band images which are coloured blue, green and red respectively. The positions of the 2MASS sources in both 
images are identified by the reference given in Table 5. 



this image has been overlaid with contours of the GLIMPSE 
8 fim and the integrated ammonia (1,1) emission coloured or- 
ange and blue respectively. This composite image reveals the 
presence of a small cluster of three very red sources towards 
the centre of Core A (these are identified in Fig. 7 by the num- 
bers 1, 2 and 3). None of these 2MASS sources are detected 
at J band, and only one is detected at H band, which would 
suggest that they are deeply embedded and that they could be 
protostellar in nature. The offsets from the centroid of the core 
position and their J, H, K s magnitudes and associated quality 
flags are presented in Table [5] 

In the right panel of Fig. [7] we present a three-colour com- 
posite image of SFO 75 produced by combining the GLIMPSE 
3.6 fim, 4.5 fim and 5.8 fim band images. Contours of the in- 
tegrated ammonia emission are overlaid in blue. The PDR can 
clearly be seen as a thin red rim of emission. This image reveals 
what appears to be an extended source which is elongated per- 
pendicular to the rim. The peak of this emission correlates with 
the brightest of the detected 2MASS sources and is extended 
in the direction of the second brightest 2MASS source (iden- 
tified as source 1 and 2 respectively in Fig. 7 and Table 5). 
The green colour is due to a 4.5 fim excess; this band is gen- 
erally dominated by stars, however, diffuse emission has been 
attributed to a combination of Brackett a and possibly H2 and 
CO in shocked regions. Diffuse 4.5 fim emission has also been 
associated with molecular outflows - it is therefore possible 
that this diffuse emission is the result of ongoing star formation 
within Core A. 



As previously mentioned, the 2MASS survey was not deep 
enough to have detected any of these sources at J band and 
only one source at H band. This limits our analysis of these 
sources and makes it impossible to use their near-IR photomet- 
ric colours to determine the evolutionary state of each source 
(i.e., protostellar, T-Tauri or ZAMS). However, we can use 
the upper limits quoted for the / band magnitude to derive a 
lower limit for spectral type and visual extinction towards each 
source. In Fig.|8]we present a K s versus J—K s magnitude-colour 
diagram. Three things are immediately apparent from this plot. 
First, two of the objects are of spectral type B2 or earlier with 
the third source being no later than B5. Secondly, all of the ob- 
jects are embedded behind ~ 18-23 magnitudes of visual ex- 
tinction. This value is approximately half that calculated from 
the H2 column density through each core (see previous sub- 
section) and places them roughly in the the centre of Core A. 
Moreover, the high values of extinction exclude the possibility 
that these sources are highly reddened main sequence stars lo- 
cated on the far side of the cloud. Finally, taking account of the 
2MASS limiting K s detection limit and assuming an average 
visual extinction of 20 mag the latest spectral t ype detectable 
would be a B5 (~ 5.9 M^: Ischmidt-Kalerll 1 9821) and therefore 
the cluster may have many more lower mass members that re- 
main undetected. 

The presence of two stars of spectral type B2 or earlier 
is in line with the findings from our mm-continuum analyses 
which suggests the presence of either a single high-mass star 
or a cluster of intermediate-mass stars. Stars with a spectral 
type of B3 or earlier are known to produce a significant pro- 
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Table 5. Parameters of red 2MASS sources located within the contoured boundary of Core A. The photometric measurements 
and their associated quality flags have been obtained from the 2MASS Point Source Catalogue (ICutri et al 1l2003l) . 



Reference 


Offset 


RA 


Dec 


2MASS Maj 


raitudes 


Quality 


ID 


(") 


(J2000) 


(J2000) 


J 


H 




flags^ 


1 


2.4 


15:55:47.3 


-54:39:11.0 


17.1 


US 


12.8 


UUB 


2 


6.6 


15:55:47.6 


-54:39:14.8 


16.8 


14.4 


12.8 


UUE 


3 


8.0 


15:55 47.1 


-54:39:16.3 


16.2 


14.6 


13.1 


UBA 



° These flags are as follows: (U) upper limit, (E) very poor fit, (B) SNR > 7 and (A) SNR > 10. 



portion of their radiation as photons with energies high enough 
to ionise the surrounding material leading to the formation of 
an ultra-compact HII region. Therefore, we might expect to 
detect radio emission associated with both of the B2 or ear- 
lier stars detected. However, no radio emission has been de- 
tected towards any of the embedded sources within Core A. 
The non-detection of any radio emission could indicate that 
the embedded stars are still at a relatively early stage in their 
evolution, somewhat later than the hot molecular core stage, 
which are generally not detected at mid-IR wavelengths, but 
before a detectable HII region has been able to form. This is 
strongly supported by the detection of a methanol maser co- 
incident with the centre of Core A recently discovered by the 
Methanol Multibeam (MMB) survey (MMB Consortium, pri- 
vate comm.fl Methanol masers are almost exclusively associ- 
ated with the earliest stages of massive star formati on turning 



off sh ortly after the formation of the UCHII region (Wa lsh et al 



1998). The detection of maser emission towards the centre of 
Core A, which has a characteristic lifetime of 2.5^1.5 x 10 4 yr 
(Ivan der Warn 120051) . allows us to place an upper limit on the 
age of the star formation taking place within Core A of no more 
than 5 x 10 4 years old. We are therefore able to conclude that 
these sources are deeply embedded massive young stellar ob- 
jects (MYSOs). 

Conversely to Core A, no 2MASS or GLIMPSE sources 
are seen towards Core B. The high values of extinction de- 
rived mean that we are limited to detecting only intermediate- 
to high-mass stars towards both cores. This makes it hard to 
rule out the presence of ongoing star formation within Core B 
completely, however, the non-detection of any stellar source 
from the 2MASS and GLIMPSE images does allow us to ex- 
clude the possibility that any massive star formation is taking 
place within Core B. If star formation is taking place within 
Core B it is likely to be limited to relatively low-mass stars 
compared to that taking place in Core A. Deeper near-IR imag- 
ing is necessary before the full stellar inventories of each core 
can be investigated in detail. 



5. Discussion 

There can be little doubt that the ionisation front is currently 
having a dramatic influence on the evolution of this cloud. 
The spatial distributions and morphologies of the IBL, PDR 
and foremost ammonia core are clearly correlated with one an- 
other, strongly suggesting these regions are interacting. In this 

8 For more details see http://www.jb.man.ac.uk/research/methanol/ 
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Fig. 8. 2MASS K s versus J—K s magnitude-colour diagram of 
the objects detected towards Core A. The thick black line shows 
K s magnitudes of main sequence stars loc ated at a distance of 
2.8 kp c, taken from lMartins & Pled d2006l) and lSchmidt-Kalei 
(119821) . and assuming zero visual extinction along the line of 
sight. The dotted lines correspond to a visual extinction of 
30 mag, w hich have bee n calculated using the standard extinc- 
tion law of MathH dl990l) with R v = 3 . 1 . 



subsection we will attempt to answer the question of whether 
the interaction between the ionisation front and the molecular 
cloud has influenced the observed star formation. 

So far we have presented detailed analyses of the structure 
of the cloud and its associated HII region. We find compelling 
evidence that SFO 75 is undergoing RDI by way of an over- 
pressured IBL and that gas dynamics indicative of shocks are 
being driven into the surface layers of the cloud. Additionally, 
we find strong evidence of recent massive star formation within 
a dense core located at the head of the cloud. The 1.2 mm 
data reveal the presence of a dense core located behind the 
cloud's rim. This core appears to be centrally condensed which 
would imply it is gravitationally bound or contains gravitation- 
ally bound objects within it, such as YSOs. The core has a 
temperature of ~ 30 K - significantly higher than found to- 
wards starless cores and Bok globules - which would suggest 
it is internally heated, further supporting the hypothesis that the 
core contains protostars that are warming the surrounding ma- 
terial. The high luminosity and mass are simila r to those found 



towards other massive star forming regions (Sridha ran et al 
20021) . 
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The higher angular resolution of the ammonia map resolves 
the 1.2 mm core into two smaller distinct components which 
are aligned perpendicular to the bright rim. It is towards the 
foremost of these two cores that we find strong evidence sup- 
porting the conclusions reached from our 1.2 mm data with 
the detection of a small deeply embedded cluster of high- and 
intermediate-mass YSOs. Since there is no evidence that star 
formation is taking place within Core B, located farther back 
from the cloud's rim, it would appear that the star formation 
is taking place exclusively very close to the edge of the rim 
of the cloud within Core A. The fact that the star formation is 
localised to the rim and is taking place within a core that is 
being strongly affected by the ionisation front and its associ- 
ated shocks is persuasive circumstantial evidence that the star 
formation has been triggered. 

The non-detection of any star formation activity towards 
Core B may indicate that while the ionisation front has had a 
dramatic affect on Core A, its associated shock front has not yet 
propagated far enough into the cloud to have had a significant 
effect on the kinematics of Core B. In this regard it is interest- 
ing to note that Core A is marginally warmer, more massive, 
and has wider line widths than Core B, in line with what would 
be expected if it is undergoing RDI. Moreover, we note that 
Core B appears to have a much simpler, almost spherical, mor- 
phology and is offset in velocity by approximately 2 km s _1 
from that of Core A. This situation is similar to examples of 
smal l scale sequential star formation seen towards other B RCs 
(e.g JSugitani et al.lll995l: lOgura et al.ll2002l:IOgurall2006l) . 

This hypothesis is supported by our analysis of the shocked 
layer identified from the CO data presented in Sect. I4.2l which 
estimated the shock had only propagated as far as the mid-point 
between the two ammonia cores. If the blue and red shifted 
emission seen in the 13 CO position-velocity diagram is due to 
the propagation of an ionisation induced shock front, we would 
expect the time taken for the shock to have propagated to this 
point to be similar to the time that the cloud has been exposed to 
UV illumination from LSS 3423. The shock propagation speed 
can be estimated using the relationship be tween the pre- an d 
post-shock pressures of the neutral gas (i.e., Whi te et al . 1999), 



V 2 (Ps ~ Pn) 

Pn 



(9) 



where P s is the pressure of the shocked gas, P n and p n are 
the pressure and density in the pre-shocked neutral gas respec- 
tively. The pressure of the shocked gas has been assumed to 
be twice that of the ionised surface layer. The value of a lies 
between 1 and 2 and using White et al's assumption that the 
pre- and post-shocked density ratio has a range between 2-oo, 
results in a maximum error of V2 into the estimated shock 
velocity. Using these assumptions we derive a shock velocity 
of ~ 2 km s , this corresponds to a propagation time for the 
shock of ~ 1.6 x 10 5 yr. 

We can estimate the amount of time the cloud has been ex- 
posed to the ionisation front by first calculating the age of the 
HII region and then comparing it to the time taken for the HII 
region to have expanded to the edge of the cloud. The bound- 
ary between the ionised gas and the surrounding molecular gas 
is clearly defined by the PDR which separates the two regions. 



We have used this interface to estimate the angular radius of 
the HII region to be ~ 3.5', corresponding to a physical radius 
of 2.8 pc. When the ionising star turns on, its radiation begins 
ionising the surrounding medium leading to the creation of a 
HII region that rapidly expands out to the Stromgren radius, 
after which point it continues to expand more slowly. Taking 
account of the Stromgren radius for a 09.5 star (0.3 pc and 
1 .5 pc for a HII region expanding into a medium of densities 
of 10 3 cm -3 and 10 2 crrT 3 respectively), and using an expan- 
sion speed of 1 1.4 km s _1 for the period of expansion after the 
Stromgren radius has been reached, we estimate the age of the 
HII region to be approximately 1-2.5 x 10 5 yr. At a projected 
distance of 0.8 pc the amount of time SFO 75 has been exposed 
to the ionisation radiation ranges between 1-2 x 10 5 yr. 

Comparing the ionisation time to the shock propagation 
time we find them to be of similar age. Although this is not con- 
clusive proof that the shock is directly connected to the arrival 
of the ionising front at the surface of the cloud, our simplis- 
tic analysis does not rule out a connection and shows them to 
be approximately coeval. Given that the star formation almost 
certainly post-dates the arrival of the ionisation front, the strong 
evidence of shocks within the surface layers of the cloud, and 
the excellent morphological correlation between the IBL, PDR 
and flattered face of Core A we conclude that the star forma- 
tion has been triggered by the photoionisation-induced shock 
driven into the cloud by the expanding HII region. The lack 
of any star formation activity within Core B further strength- 
ens the case for the star formation within Core A to have been 
triggered. 

6. Summary and conclusions 

We present a multi-wavelength study of the bright-rimmed 
cloud SFO 75 in order to evaluate whether or not the observed 
star formation could have been triggered by radiative-driven 
implosion. This study includes observations of the 1.3 cm and 
1.2 mm continuum, and I3 CO (7=1-0) and ammonia (1,1) 
spectral line observations. To complement our observations 
we have used near-, mid- and far-IR archival data obtained 
from the 2MASS, GLIMPSE, and IRAS archives respectively. 
Combining all of these data sets has enabled us to build up a 
comprehensive picture of this cloud and its local environment. 

RCW 98 appears to be an approximately spherical HII re- 
gion surrounded by a low density (< 10 3 cm 4 ) ring of molec- 
ular material. The HII region has an angular radius of ~ 3'. 5, 
corresponding to a physical radius of ~ 2.8 pc, from which we 
estimate its age to be 1-2.5 x 10 5 yr. SFO 75 is located on the 
south-east side of the HII region, and approximately lies in the 
same plane in the sky as the ionising star. The bright rim of 
SFO 75 is located ~ 0.8 pc from the ionising star and has been 
exposed to the ionising radiation for ~ 1-2 x 10 5 yr. 

We are able to clearly identify the ionised boundary layer 
(IBL), the photo-dominated region and the molecular material. 
The excellent morphological correlation between these three 
leads us to conclude that there is a high degree of interaction 
between them. The IBL is over-pressured with respect to the 
molecular gas by a factor of three and it is therefore likely that 
photoionisation-induced shocks are currently being driven into 
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this cloud. This is strongly supported by the detection of red 
and blue shifted emission which is attributed to emission from 
shocked gas. We estimated the shock propagation time to be 
~ 1.6 x 10 5 yr and find it to be approximately coeval with the 
cloud's exposure to the HII region. 

Our CO and 1.2 mm data reveal the presence of a dense 
molecular core located slightly behind the bright rim. The den- 
sity and mass of the molecular gas range between ~ 0.3- 
3 x 10 5 crrT 3 and 170-570 M respectively. We find strong 
evidence of recent massive star formation within a dense core 
located at the head of the cloud. This core appears to be cen- 
trally condensed which would imply it is gravitationally bound 
or contains a gravitationally bound object within it. The core 
temperature (~ 30 K) is significantly higher than found towards 
starless cores and Bok globules indicating the presence of an 
internal source of heating, such as YSOs. Moreover, the high 
luminosity and mass are simi lar to those found tow ards other 
massive star forming regions ( Sri dharan et al ] l2002h . 

The 1.2 mm core is resolved by our high resolution am- 
monia observations into two distinct cores, Core A which is 
located directly behind the bright rim and Core B which is lo- 
cated farther back. Towards Core A we find a small cluster of 
MYSOs (two B2 and one B5 stars or earlier) which are deeply 
embedded behind ~ 20 mag of visual extinction. The upper 
limit for the age of this cluster, which is approximately half 
the estimated ionisation time, would indicate that its formation 
post-dates the cloud's exposure to the HII region. Due to the 
limited magnitude of the 2MASS observations we are unable 
to detect any sources that might be present that are later than 
B5, so these sources might just be the most massive members 
of a larger cluster that remains hidden. No stars are found to- 
wards Core B and given its simple morphology and that it is 
located farther back from the rim than Core A leads us to con- 
clude that the shock front has not yet penetrated far enough into 
the cloud to have reached this core. 

Although we are unable to determine conclusively whether 
the observed star formation has been triggered we have pre- 
sented a considerable amount of circumstantial evidence which 
we believe strongly supports the hypothesis that the star for- 
mation has been induced. The scenario that emerges from our 
analysis is one where the two ammonia cores pre-date the ar- 
rival of the ionisation front. However, since its arrival the over- 
pressure of the ionised gas at the surface of the cloud has driven 
shocks into the surface layers of the cloud. The propagation of 
these shocks through Core A has triggered the formation of a 
small cluster of massive stars, however, the shock front has not 
yet propagated deeply enough into the cloud to have affected 
the evolution of Core B. 
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